Quantum phase diagrams of fermionic dipolar gases for an arbitrary orientation of 

dipole moment in a planar array of ID tubes 
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We systematically study ground state properties of fermionic dipolar gases in a planar array of 
one-dimensional potential tubes for an arbitrary orientation of dipole moments. Using the Lut- 
tinger liquid theory with the generalized Bogoliubov transformation, we calculate the elementary 
excitations and the Luttinger scaling exponents for various relevant quantum orders. The com- 
plete quantum phase diagrams for arbitrary polar angle of the dipole moment is obtained, including 
charge density wave, p-wave superfluid, inter-tube gauge-phase density wave, and inter-tube s-wave 
superfluid, where the last two breaks the U{\) gauge symmetry of the system (conservation of 
particle number in each tube) and occurs only when the inter-tube interaction is larger than the 
intra-tube interaction. We then discuss the physical properties of these many-body phases and their 
relationship with some solid state systems. 
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I. INTRODUCTION 

It is well-known that the system of ultracold atoms has 
been the most experimentally flexible system to study the 
strongly correlated phenomenons in the recent decade, 
due to the widely tunable interaction strength, dimen- 
sionality and/or the species of underlying atoms. Since 
the many-body physics in one-dimensional systems has 
been extensively studied by various analytic and numeri- 
cal methods in the last few decades, it is therefore impor- 
tant to compare the theoretical predictions to the exper- 
imental observation in the recent low-dimensional cold 
atom systems. Several important works along this direc- 
tion have been reported recently, including the Tonks- 
Girardeau gas fH , Luttinger liquid (LL) behavior [2] , two- 
component fermionic gas 3], ID BEC-BCS crossover 
polaronic effects in Bose-Fermi mixture Q, and ID spinor 
gases @ etc.. However, due to the short-range nature of 
atomic interaction, it is usually not easy to study how the 
interaction between particles in different ID tubes can 
bring different many-body effects, which, however, can 
be easily achieved in the traditional solid state system 
due to the long-ranged nature of Coulomb interaction. 
Knowing the fact that the dipolar interaction between 
spinful atoms and polar molecules can be also prepared 
in various circumstances 0@| , it is therefore getting more 
and more attention on the effect of long-ranged dipolar 
interaction between ultracold atoms/molecules in low- 
dimensional systems [olfioj. Among these earlier works, 
we are especially interested in the systems of multi-tubes 
P ITllfl^ . which provides a direct example to study how 
the quasi-lD system can be merge to the physics of two- 
dimensional system. 

In this paper we extend our earlier work in the study of 
double-tube systems [llj to investigate the ground state 
properties of a planar an array of ID tubes loaded with 
dipolar fermionic atoms/molecular, assuming the inter- 
tube tunneling is energetically negligible. The dipole mo- 
ments are polarized by an external electric or magnetic 
field in an arbitrary direction (see Fig. [1]), providing a 




FIG. 1: (Color online) A planar array of ID tubes of a dipolar 
gas with the inter-tube distance, d, the tube radius, R, and 
the tube length of tube, L. In this paper, we assume L S> 
R, and therefore the system can be regarded as a quasi- 
lD system. The dipole orientation is along the direction of 
the external field, E, which is tilted with a polar angle 9 from 
the z axis. The azimuthal angle is denoted by (p. 



variety of interaction matrix elements that are not avail- 
able in traditional solid state systems. Using the multi- 
component Luttinger liquid theory, we then diagonalize 
the effective system Hamiltonian exactly, and obtain the 
quasi-long ranged quantum orders by studying the cor- 
relation functions. We further study the elementary ex- 
citations and investigate the regime of system instability 
via a mode softening in the long wave-length limit. 

Our results show that, some unexpected inter-tube 
correlation or pairing mechanism can be found when 
the dipole moment is tilted near a magic angle, 9c = 
cos^^ (see Fig. [T|), where the intra-tube interac- 

tion is relatively smaller than the inter-tube interaction. 
Depending the inter-tube interaction is positive or neg- 
ative (for different angle, (/), see Fig. [1]), fermions in 
the neighboring tubes can have inter-tube gauge-phase 
density wave (GPDW) or inter-tube s-wave superfluid 
(s-SF), breaking the U{1) gauge symmetry (i.e. par- 
ticle number conservation in each tube) of the system 
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(in the quasi-long-ranged order sense). We note that in 
the quantum Hall double layer system, the gauge sym- 
metry breaking phase (or called inter-layer coherence, 
pseudo-spin ferromagnetism, or exciton condensate) has 
been observed experimentally in a semi-conductor het- 
erostructure , showing an interaction-enhanced single 
particle zero-bias tunneling current. We believe similar 
results may be also observed in the multi-tube system we 
discuss here. For the inter-tube s-wave superfluid state, 
we expect a perfect transmission for a particle moving 
through a potential barrier along the longitudinal direc- 
tion of the tubes at zero temperature, and a power-law 
decaying behavior at finite temperature, similar to the 
results predicted by Kane and Fisher [ij] on a single in- 
teracting ID electronic gas. When the angle 9 is larger 
or smaller than 9c, intra-tube interaction becomes domi- 
nant and the system becomes simple charge density wave 
(CDW) or intra-tube p-wave superfluid (p-SF) phases, 
consistent with the standard mean field results. (We will 
show a clear definition and physical picture for each of 
these phases later.) Our works therefore provide a clear 
evidence of the strong dipolar long-ranged interaction ef- 
fects in this quasi- ID systems. 

This paper is organized as following: In Sec. |TT1 
we present the interaction matrix element, the system 
Hamiltonian in Luttinger liquid theory, and the diagonal- 
ization method for the multi-component Luttinger liquid 
Hamiltonian. Wc then study the elementary excitations 
in Sec. Illll and the correlation functions in Sec. IIVI where 
some details of calculation are shown in Appendix [X] Wc 
then show the complete quantum phase diagram for arbi- 
trary dipole angle and discuss related physical properties 
in Sec. IVl We then summarize our results in Sec. IVII 



II. SYSTEM HAMILTONIAN AND 
DIAGONALIZATION METHOD 

A. Interaction matrix element 

The multi-tube system considered in this paper is 
shown in Fig. [H where the trapping potential can be 
generated by a 2D optical lattice in the x — y plane with 
a stronger lattice potential in the x-axis. Additional weak 
magneto-optical potential is applied in the longitudinal 
direction (z) to trap dipolar atoms/molecules in such 
quasi-lD potential. In this paper, we assume the lon- 
gitudinal length scale, L, is much larger than the inter- 
tube distance d and tube radius R (i.e. L ^ d ^ R). 
For simplicity, we assume the inhomogeneity of the weak 
trapping potential along the y and z axises can be ne- 
glected, and therefore the system can be regarded as an 
array of identical ID tubes with total number of tubes, 
N. 

Throughout this paper, we always consider the case 
when dipolar atoms/molecules are all loaded in the lowest 
subband of each tube, which has a transverse confinement 
wavefunction: 0j(x,y) = -A_e-to'+(a-i<i/2)2)/2fl2_ -pj^g 



resulting effective ID system Hamiltonian then can be 
written to he H = Hq + Hj, where Hq is the kinetic 
energy: 



^^0 = 



E 



(1) 



with m being the mass of dipolar particles and tpjiz) 
being the fermion field operator in the jth tube, and Hj 
is the interaction energy: 
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Here V\j^j>\{z) is the dipolar interaction between 
molecules in the same tubes (j — j') or in different tubes 
[i ^ j''), obtained by integrating out the transverse de- 
gree of freedom (rj^ = (x,y)): 

V\j_j.\{z) = j dria j dr2,_L |(/)j(ri,_L)|^ |(Aj'(r2,j_)|^ 

xt/d(ri-r2), (3) 

where Vd{v) = D'^{\ - 3(f • Ef )/\v\^ is the bare dipo- 
lar interaction with D being the electric dipole moment 
in the c.g.s unit. E is the unit vector along the exter- 
nal field, parallel to the direction of dipole moment (see 
Fig. [T|). Since in general the electric dipole interaction 
is much stronger (and tunable) than the magnetic dipole 
interaction, in the the rest of this paper we will use polar 
molecules as the underlying particles for further discus- 
sion. Extension of our results to the magnetic dipolar 
atoms is straightforward. 

For the convenience of later study, here we define 
the interaction matrix element to be the Fourier trans- 
form of the bare dipole interaction in momentum space, 
i.e. we take the first Born approximation for the two- 
particle scattering for simplicity (see Ref. [IBl): we have 
V\j_ji\{q) = / dzV|j_jv| (z) e^*'^^. Although one has to 
evaluate Mj-j'|('?) numerically for finite value of q, we 



can still obtain an analytic form of their zero momentum 
{q = 0) values, which dominates the low energy physics 
in the Luttinger liquid theory. By integrating over the 
transverse confinement wavefunction, we have 
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for the intra-tube interaction and 
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for the inter-tube interaction (I = 
to see that when 9 = 9c = cos" 



b'-j'D- It is very easy 
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FIG. 2: (Color online) The inter-tube interaction for two 
neighboring tubes, \j — j'\ — 1, with arbitrary dipole ori- 
entation and d/R = 5, see Eq. ([5}. 



FIG. 3: (Color online) Different scattering process for 
Tomonaga-Luttinger model. We assume the system is away 
from half filling regime and consider the small momentum 
transfer only, so that g^^^ and g^'^^ scatterings are neglected 
in our system. Therefore, we only consider forward scattering, 
g'^-* and g''*', in our system Hamiltonian. 



interaction becomes zero and change sign, while the inter- 
tube interaction vanishes and changes sign a,t ip — tt/A. 
In other words, we can have a wide parameter range to 
tune the sign of interaction matrix element by tuning 
the direction of external field to explore various kinds 
of interesting many-body physics. In Fig. [21 we show 
the value of nearest neighboring inter-tube interaction 
matrix element, Vi(0), as a function of 9 and (j). We 
note that such interesting kind of inter-tube interaction 
matrix element cannot be realized in the traditional solid 
state system and therefore may bring some physics even 
not predicted before in the condensed matter theory. 



B. Luttinger liquid Hamiltonian 

For one-dimensional non-interacting fermions, the ele- 
mentary particle-hole excitations around the two Fermi 
points can be well-approximated by the linearized disper- 
sion, known as the Tomonaga-Luttinger model [mill. 
In the low energy limit of an interacting system (i.e. 
when the characteristic energy scale is much smaller than 
the Fermi energy), the most important excitations are 
still around the Fermi points and therefore the ground 
state and the elementary excitation properties can be 
well-described by the so-called Luttinger liquid theory 
' T8lfl9ll20| ■ where the single particle excitation of the 



Landau-Fermi liquid becomes totally destroyed by the 
inter-particle interaction, leaving only gapless collective 
modes as the low energy elementary excitations. Within 
the standard Luttinger liquid theory, only four kinds of 
scattering near the two Fermi points are important (see 
Fig. O: .9 (2) and (7^^) are the forward scattering without 
changing the direction of scattered fermions, g^^^ is the 



backward scattering with the change of momentum, 2kp, 
and g^'^^ is the Umkalpp scattering, which is important 
only when a lattice potential is present and the filling 
fraction is close to unit. In our present work, we are 
interesting in the regime when g^^-* is negligible or irrel- 
evant, and (7^^^ can be omitted since no lattice potential 
along the longitudinal direction of ID tubes is consid- 
ered here. Renormalization group for these multi-tube 
systems can be used to study the relevance of the back- 
ward scattering term (i.e. g^^-*), and the results when it 
becomes relevant will be presented in another place in 
the future plj. Throughout this paper, we are interested 
in the regime when this backward scattering is irrelevant 
or at least negligible. 

In the multi-tube system we consider here, the two 
scattered particles can be either in the same tube (intra- 
tube interaction) or in two different tubes (inter-tube in- 
teraction). As a result, labeling the single particle quan- 
tum states by the tube index, j {j — 1, - ■ ■ , N), and the 
left /right mover index near the two Fermi points (r = ±), 
the system Hamiltonian of Eqs. ^ and ^ can be writ- 
ten as 2N component Luttinger liquid model: 



Hll = 



E 

id' 



H 



(2) 



H 



(4) 



(6) 



where j, j' are the tube index. The noninteracting Hamil- 
tonian, -ffj"^ , can be approximated by a linearized band 
structure about the Fermi points: 

r— ± q 

9>0,r=± 
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Here Cj^r,q is the Fermionic creation operator for dipoles 
with momentum q. vp and kp are the Fermi vefocity 
and Fermi momentum, and are assumed to be the same 
for all tubes. In the last line, we have used the standard 
result of Luttinger model, where the kinetic energy of 
a linearized dispersion can be composed of two density 
operators [2^, pj,r{q), where 

PjAq) = 4:r,fc+gCj,r,fc. (8) 
k 

Note this density operator here obeys a boson-like com- 
mutation relation, 

= 5q,q'6x.y^sign{vx), (9) 

where we define A — {j, r) to label the tube and char- 
ilty indecies, and v\ = rvp. As a result, v\ = :kvp for 
the right/left-moving channel. Such convention can be 
very useful for our latter application to the generalized 
Bogoliubov transformation. 

As for the interaction Hamiltonian of forward scatter- 
ing, it is easy to show that it can be also rewritten to be 
the following form within the Luttinger liquid theory: 

= \^ll9f}M)p,Aq)prA~<D (10) 
- ^ E af}^ii)pU<DPrA-<i)- (11) 

q,r—±. 

Note that in the limit of weak interaction (i.e. interaction 
energy is smaller than the Fermi energy) the two kinds 
of forward scattering (g^^^ and 5*^^-') are the same as the 
bare interaction, i.e. 

9^'^ = 9f}'{<l) = V\,-rM)- (12) 

Since both the noninteracting Hamiltonian and the inter- 
action Hamiltonian are now quadratic in terms of density 
fluctuation operator, the total Luttinger liquid Hamilto- 
nian, Eqs. ([7]), (fTO|) and (fTT|) . can be diagonalized by us- 
ing a generalized Bogoliubov transformation as described 
below. 



C. The Generalized Bogoliubov transformation 

To make the notation more compact for the Bogoliubov 
transformation of a general A'^-tube system, we use A to 
denote both the tube index (j) and the chirality index 
(r — ±). The whole Hamiltonian shown above can be 
described to be 

^=Y E Ax,w[q)px{q)px'[-q). (13) 



where 

Ax.x'{q) = \vx\5x.x' + ^^^^^ (14) 

ZTT 

is the matrix element of the matrix A. 

We diagonalized Eq. ([13]) via the generalized Bogoli- 
ubov transformation developed by Penc and Solyom [2^ . 
Defining the eigenstate density fluctuation to be Pn{(i) 
with the index rt = 1, 2, • • • , 2A^, the transformation can 
be written to be 

Pn{q) = X! ^»,aPa(9), (15) 

A 

where Wn,x is the matrix element of a canonical transfor- 
mation (i.e. the generalized Bogoliubov transformation) 
matrix. After the transformation, the Hamiltonian in Eq. 
(fT3)) should be expressed to be 

^ = Y"EE l""l^"(^)^"(~^) 

q>0 n 

= X! Y^n^\Px{q)\Un\Wn^X'P\'i~q) (16) 

i3>0,?i A, A' 

In other words, comparing Eq. (fT^ and Eq. (fTB)) . we 
must have 

^A,A' = ^ \Un\Wn,XWn,X' (17) 
n 

or 

A = ^|u„||ii;("))(z«(")|, (18) 

n 

where jw^"^) is a vector with elements, Wn,x- The com- 
mutation relation of the new density operators can be 
also expressed to be (following Eq. ([9])): 

[pn{-q),Pn'{q')] = <5,,,.^(u;(")|B|u;("')), (19) 

ZTT 

where the matrix B in these two different basis should 
be: 

Bx^X' - Sx^x'Sign{vx) (20) 
B„,„. EE = ,5„,„'sign(u„). (21) 

From the above derivation, we want to find a basis 
{jw*^"-')} to diagonalize the matrix A with positive eigen- 
values, |u„|, and the eigenvectors have to satisfy the spe- 
cial orthogonal relation in Eq. ([2T|) . It is easy to see that 
it is equivalent to solve the standard eigenvalue problem 
below: 

AB|w(")) =w„|w(")) (22) 

with the condition of Eq. pT|) . This is the so-called 
generalized Bogoliubov transformation for the multi- 
component Luttinger liquid system, first developed by 
Penc and Solyom 
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q//cp q//cp 

FIG. 4: The calculated energy dispersion for the eigen-mode 
excitations of four tubes system, (a)-(d) are for 7 = 2, 6, 10, 
and 14 respectively, and we set kpd — 7, dn = 5, 9 = 6c, and 
= here. The instability occurs at the lowest branch for 
7 = 14. 



III. ELEMENTARY EXCITATIONS AND 
SYSTEM INSTABILITY 

A. Results for = 

The first physical properties we want to study is the el- 
ementary excitations of such multi-tube system. For the 
convenience of later discussion, here we introduce sev- 
eral dimensionless parameters to describe the systems: 
7 = is to measure the dipolar interaction strength, 
dj^ = d/R is the inter-tube distance compared to the tube 
radius, and kpd is to measure the particle density in each 
tube. Since the long wave-length intra-tub e in teraction 
can be tuned to zero as ~ 6*0 = cos~^ a/1/3 (see Eq. 
(|3])), it is reasonable to investigate the results for 9 = 9c 
first. We also set = 0, i.e. the dipole moment is in the 
x—z plane first for the convenience of study. Results for a 
more general angular parameter regime will be presented 
later. In Fig. 21 we show the calculated dispersion of col- 
lective modes for different values of interaction strength, 
7, in the system of four tubes (iV = 4). It is easy to see 
that since the intra-tube interaction is zero at 9 = 9c, the 
dispersion of the collective excitation becomes hybridized 
by the inter-tube interaction, which lifts the degeneracy 
and changes the excitation dispersion dramatically in the 
long wavelength limit, while it becomes less important in 
the short wavelength (large momentum) regime. The up- 
per two branches exhibit convex shape at small momen- 
tum and the lower two branches show a concave shape. 
This train is enhanced in strong dipole-dipole interaction 
regime, and the lowest energy branch eventually becomes 
softened when the dipole strength, 7 is above a critical 
value. By investigating the eigenstate wavefunction of 
the lowest energy collective mode, we find that the insta- 
bility should be related to some kinds of phase separation 
state, which is derived by the strong inter-tube repulsive 
interaction in such multi-component system. 

In Fig. [Ua)-(c), we show how the velocity of collective 




FIG. 5: (Color online) (a)-(c) show the group velocity at q = 
as a function of interacting strength 7 for A'^ = 2, 3, and 
4 respectively. Here we set dn = 5, 9 — 9c, and kpd = 
7. (d) shows the critical value, 7c, of the softening of the 
lowest energy excitation as a function of numbers of tubes, 
A'^. Results for kpd = 1, 2, ■ ■ ■ 7 are shown together and found 
scaled to a single universal curve, above which the system 
becomes unstable. 




q/kp q/kF 

FIG. 6: Elementary excitation energy dispersion for 9 — 9c, 
cj> — 7r/4, and kpd = 7. 7 = 2, 6, 10, and 14 respectively from 
(a) to (d). 

modes, Uj, changes as a function of interaction strength 
for different numbers of tubes. Defining jc to be the 
critical interaction strength for the softening of the lowest 
excitation energy, in Fig. El^d) we show that this critical 
value is proportional to kpd upto a universal function of 
the number of tubes, N. The physical interpretation of 
such result is straightforward: the ratio between 7 and 
kpd is equivalent to the ratio of interaction energy to 
the Fermi energy, which is the only important parameter 
to determine the quantum phase diagram in our present 
system since the intra-tube interaction is zero here. 

B. Results for / 

When 7^ 0, the dispersion relation will change in an- 
other fashion. In Fig. [51 we show the energy dispersion at 
= 7r/4. At this angle, both inter- and intra-tube inter- 
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jt/2'O 



FIG. 7: (Color online)The critical interaction strength for the 
system stability as a function of dipole moment polar angle, 
9 and 0. The number of tubes is A'^ = 8. We show the results 
both in a 3D plot and the 2D shading plot. The red and blue 
lines respectively the line with 9 = 6c and — 7r/4. 



action matrix elements are zero in the long wave-length 
limit, i.e. ^/(O) = for Z 0, 1, 2 • • • . The systems be- 
comes different from noninteracting systems only through 
the interaction at finite momentum regime {k ^ 0). We 
can see that when the interaction 7 changes from from 
2 to 14, the strong interaction separates the dispersion 
curves in a totally different way from (j) — case. Dif- 
ferent branches of excitations intersect with each other 
at none zero momentum, and the intersection point is al- 
most the same even when the dipole strength increases. 
When the interaction is stronger than a critical value, 
the system becomes unstable through a softening at fi- 
nite k as expected. In Fig. [71 we show the calculated 
critical interaction strength for = 8, as a function of 
6 and 0. We can see that the system becomes more sta- 
ble when 9 is larger than the magic angle, 9c, while the 
variation of angle (p does not affect the stability of the 
system very much, except for some asymmetric critical 7 
between (f> > tt/A and (p < tt/A. 



IV. CORRELATION FUNCTIONS AND 
SCALING EXPONENTS 

Now we study the correlation functions within the Lut- 
tinger liquid theory. The correlation functions of the or- 
der parameters O should decay in a power-law in the 
large distance limit [l9ll25ll26j : 



(on^)o(o)) 



1 



(23) 



where a is the associate scaling exponent. Here we have 
assumed that the order parameter, O, here is composed 
by two-fermion operators only and therefore an order pa- 
rameter is of quasi-long-ranged order only when a > 0. 
The order parameters of the largest Luttinger component 
should be understood as the dominant quantum phase 
for a given parameter. The detail calculation of the Lut- 
tinger exponent is shown in Appendix. [X] 

Based on our earlier results in the double tube system 
of Ref. here we can just consider the following four 
kinds of order parameters: charge density wave (CDW), 
p-wave superfluid (p-SF), inter-tube gauge phase density 
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FIG. 8: Luttinger exponents for different order parameters 
a,t (j) — 0. for the GPDW phase, we calculate the correlation 
function between the first tube and the other tubes, j, labeled 
by the horizontal index. For the other two phases, we show 
result of correlation function of a single tube. Here we set 
7 = 2, dfl = 5, kpd = 7, and iV = 6. 9 = 9c - 27r/36, 
9c ~n/36, 9c, Sc-l-7r/36, and 6c + '2ti respectively from (a) 
to (e). 
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FIG. 9: Luttinger exponents for different order parameters 
at = 7r/2. for the s-SF phase, we calculate the correlation 
function between the first tube and the other tubes, j, labeled 
by the horizontal index. For the other two tubes, we show re- 
sult of correlation function of a single tube. Other parameters 
are the same as Fig. [H] here 6 — 6c ~ 2n/36, 9c — tt/36, 6c, 
9c + vr/36, and 9c + 27r/36 respectively from (a) to (e). 
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wave (GPDW), and inter-tube s-wave superfluid (s-SF). 
Note that we redefine the names of some order parame- 
ters, because the usual Luttinger liquid theory has only 
two (pseudo-spin) components, while the case we stud- 
ied here has more components. The (pseudo-)spin based 
terminologies (say, triplet/singlet pairing or spin density 
wave etc.) are therefore not suitable here at all. How- 
ever, we can still define the order parameters in any pair 
of tubes as following: 

O^tsFix) ^ 4_(x)4+(a;) (26) 
Ol_sA^) ^ 4+(x)4__(a;)-4^+(x)4_(x)(27) 

We also have used the fact that the dominant quantum 
phases are always contributed from fermions of two op- 
posite chiralities (r = ±). 

In Fig. [H we show the calculated Luttinger scaling ex- 
ponents for each of these order parameters at = 0. 
Since s-wave superfluid never becomes a candidate of 
quasi-long-ranged order (i.e. its exponent, a, is always 
negative in the entire range we consider here), we will 
not show its result for simplicity. Note that there is no 
correlation between fermions in tube j and tube j' for the 
CDW and p-SF phase, and therefore their Luttinger ex- 
ponent, a, is finite only when j = j' . As a result, in Fig. 
|8]we show the obtained exponent for each tube (j — j') 
if considering these two phases, while we show results 
between tube j = 1 and j' ^ 1 when considering the 
GPDW phase (similarly for the s-wave superfluid phase, 
but not shown here). We can find that for 9 < 6,., the 
dominating phase is p-wave superfiuid, due to the strong 
attractive interaction between fermions in the same tube. 
Increasing 9 gradually shows that the this phase becomes 
suppressed, and the inter-tube gauge phase density wave 
(GPDW) becomes dominate when 9 ~ 9c. For 9 is get- 
ting larger, the ground state becomes dominated by the 
charge density wave (CDW) phase eventually. In a dou- 
ble tube system (A^ = 2, Ref. [lT'|), such GPDW phase is 
also called the planar (pseudo-)spin-density wave phase, 
breaking the gauge symmetry of the particle conservation 
in each tube. In the mean field level, this corresponds 
to an interaction induced effective tunneling correlation 
between neighboring tubes, similar to the ferromagnetic 
state in the quantum Hall bilayer systems p^, except 
that the gauge phase here is not uniform, but oscillating 
along the tube. 

In Fig. [HI we show results in another limit, = 7r/2, 
which makes the inter-tube interaction attractive. Again, 
we can see that the p-SF and GDW phases are the domi- 
nant phases for 9 is smaller and larger than 9c. However, 
when 9 ^ 9c, the s-wave superfluid becomes dominant, 
showing an inter-tube pairing phase. Different from the 
typical BCS pairing, here the pairing is between any two 
nearest neighboring tubes and therefore can be regarded 
as a special multi-component BCS pairing phase. 
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FIG. 10: (Coloer online) Typical quantum phase diagram for 
multi-tube systems for kpd = 7, dn = 5 and 7 — 0.6. Right 
hand side are cartoons for the four dominant phases in the 
system. The black filled circles indicate fermions and the open 
circles indicate holes. The elliptic circles with solid/dashed 
lines indicate pairing/coherence between the two fermions. 



V. COMPLETE PHASE DIAGRAM AND 
DISCUSSION 

Combining above results we can summarize the ground 
state properties of weakly interacting dipolar fermions in 
a system of a planar array of ID tubes in Fig. [TOl we 
also show some cartoon pictures to describe the physical 
meaning of each relevant order parameters in this sys- 
tem. We can see that by tuning the polar angle of dipole 
moment, we can have these four different phases (CDW, 
GPDW, p-SF, and s-SF) in different regime of the phase 
diagram. The main features to make this system different 
from traditional solid state systems are the appearance 
of inter-tube gauge-phase density wave (GPDW) and the 
s-wave superfluid (s-SF), where the former requires a 
stronger inter-tube repulsive interaction compared to the 
intra-tube interaction, while the later requires a stronger 
negative inter-tube interaction. Both of these two phases 
break the C/(l) gauge symmetry of this system (conserva- 
tion of particles in each tube) . Since the CDW phase and 
the intra-tube p-wave superfluid phase are easily under- 
stood from the nature of long-ranged repulsive and at- 
tractive interaction, here we concentrate on the physical 
picture and implication of the other two phases, inter- 
tube gauge phase density wave (GPDW) and inter-tube 
s-wave superfiuid (s-SF) phases. 

The existence of the GPDW phase has been observed 
in a double-tube system in Ref. where the pseudo- 
spin language is used to define it as a planar pseudo- 
spin density wave. As have been mentioned in Ref. 
such GPDW phase breaks the f/(l) symmetry of the un- 
derlying Hamiltonian (i.e. particle conservation in each 
tube), similar to the quantum Hall ferromagnetism in a 
double-layer system [13|. In the quantum Hall system, 
the pseudo-spin ferromagnetism (also called inter-layer 
coherence or exciton condensate) is experimentally man- 
ifested by an interaction induced tunneling current in the 
limit of zero single particle tunneling. In our present 
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system, we have such ID inter-tube correlation between 
any two nearest neighboring tubes, and therefore beheve 
that if a small (but energetically negligible) single parti- 
cle tunneling is allowed, the tunneling current (quantum 
fluctuations) between any two nearest pair of tubes can 
be greatly enhanced, showing a nontrivial long-ranged 
correlation in the direction perpendicular to the tube di- 
rection also (see the double-tube result in Ref. [IH)- 
In other words, a dipolar atom or molecule in the first 
tube has a finite probability to be found in the tube of 
the other side, even the single particle tunneling rate is 
much smaller than the value needed for the same result 
in the noninteracting case [28j . Such interesting many- 
body effect may be also relevant to the Luttinger liquid 
explanation for the high- Tc superconductor in the two- 
dimensional cuprate, proposed by P. W. Anderson [29(. 
We believe it is worthy to investigate further along this 
direction by using systems of ultracold polar molecules 
as a quantum simulator, which can have a very strong 
dipolar interaction to make the scaling exponent, a, to 
be large enough to be observed experimentally. 

As for the inter-tube s-wave superfluid, we empha- 
size that it is different from the usual BCS type pair- 
ing by additional multi-component physics: the Cooper 
pairing can occur between any two nearest neighboring 
tubes. For a two-component ID system (say, spin-half 
electron gas), the attractive interaction between fermions 
also leads to the s-wave pairing and hence a perfect tun- 
neling through a potential barrier at zero temperature, 
as first investigated by Kane and Fisher [13] by using 
renormalization group method. At a finite temperature, 
they also predicted a power-law decrease of the tunneling 
rate. Therefore, it is then reasonable to expect similar 
behavior should be also observed for the pairing between 
multi-component fermions. Moreover, since there are at 
least two choices for a fermionic dipole in the middle tube 
to pair another dipoles in its two neighboring tubes, such 
superflow property should be more significant than the 
result of two-component pairing case. A more extensive 
study of the transport property of such multi-tube sys- 
tems can be very interesting and significant both in the 
theoretical and experimental sides, but has beyond the 
scope of this present paper. 



VI. SUMMARY 

In summary, we study the complete quantum phase 
diagram of a fermionic dipolar atoms / molecules in a pla- 
nar array of one-dimensional tubes. Using the Luttinger 
liquid theory and the generalized Bogoliubov transforma- 
tion, we are able to derive the elementary excitations and 
the Luttinger components for various correlation func- 
tions. From the mode softening of the lowest energy ex- 
citations, we find that such a fermionic dipolar system 
will become unstable when the interaction is stronger and 
when the number of tubes increases. From the calcula- 
tion of Luttinger exponents, we identify the parameter 



regimes for various kinds of order parameters to be dom- 
inant, including charge density wave, p-wave superfluid, 
gauge phase density wave, and s-wave superfluid, where 
the last two exist only when the intra-tube interaction is 
smaller than the inter-tube interaction. These two inter- 
esting many-body phases cannot be realized in the semi- 
conductor based quantum wire systems, and are worthy 
to be investigated further in the future theoretical and 
experimental studies. 

The authors appreciate valuable discussion with H.-H. 
Lin, C.-M. Chang, M. Cazalilla, and E. Demler. This 
work is supported by NSC (Taiwan). 



APPENDIX A: LUTTINGER SCALING 
EXPONENTS OF CORRELATION FUNCTIONS 

In this section, we will briefly introduce the method to 
calculate the scaling exponents of various order param- 
eters within the bosonization method. We first express 
the fcrmion operator by the following bosonized expres- 
sion t2Qi26|l: 



lim —^= 

"^0 y'2TTa 



exp [-i{r^x{x)-ex{x))] (AI) 



where a is a short-ranged cut-off and will be taken to 
be zero in the final expression of any physical quantity. 
We also have neglected the fermionic creation operators, 
which is not relevant to the calculation of correlation 
function here. The two new bosonic fields are defined 
to be 



E 

q^Q 



p-0!\q\/2^~iqx 



(A2) 



ITT X ^ 

q^O 



-a\q\/2 -iqx 



[Pj,+ i<l) - P3,-i<l)] : 



(A3) 



which are the bosonized density fluctuation and phase 
fluctuation operators respectively. In this section, we will 
use j to label the different tubes and therefore A — {j, ±) 
is to label the fluctuations of the jth tube with chirality 
r = ± (i.e. the right/left movers). To transform the bare 
density fluctuation to the eigenstate bases, pn{q), we use 
the following inverse transformation for Eq. ([TS]): 



2N 



Pj 



11=1 



2N 



Pj, 



-(<7) = ^Tj+iy^„pn{q), 



(A4) 



(A5) 
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where n is the eigenvector index, and n < iV for right 
moving mode {r — +) and n > N+1 for left moving mode 
(r ~ — ). For symmetric model, i.e. the Fermi velocities 
of the left movers and the right movers are the same, the 
matrix, T, will be block symmetric, i.e. Tj^n — Tj^j^_n+N 



and T, 



j+N,n 



for n < A^. We can insert Eq. (|A4|) 



and Eq. ((X5|l into Eqs. (A2^ and ((X3|l and simphfy the 
result using the block symmetry. We obtain 



N 



ITT f.-a\l\/-2^-rqx 



[Pnil) +Pn+Niq)] 



N 



N 



{Tj,n + T,^N,n) ^n{x) EE ^ C^„|.„(x) (A6) 



n=l 



where 



Fi{x) = (($„(a;) - $„(0))2) = ((e„(x) - e„(0))2XAl2) 



as defined in the Ref. [26|. Note that, we have used 
the fact that the CDW and p-SF phases are defined for 
fermions in the same tube (j — j'). The universal func- 
tion, Fi{x), has been calculated in details in Ref. [2^ . 
and here we just show the final result: 



■In 



(A13) 



N 



Y E : [^"('^) ^ Pn+N{q)] 



q 



N 



N 



(A7) 



n=l 



where = T,,„ + Tj+N,n and = Tj^n - Tj+N,n- 

The next step is to use the bosonized fermion field 
operator in Eq. (jAl[) to calculate the desired correla- 
tion function, Eqs. (|24 p - (|77|) . We can combine all the 
phase fields into one exponential by using the following 
identities: e^e^ = e'4+-Bg[A,B]/2^ ^^^^ ^^A^ ^ g^i^') for 
bosonic operators. After some straightforward calcula- 
tion, we find that all the correlation functions we con- 
sider in this paper can be simplified into the following 
forms 129: 



oc exp 



(A8) 



i^G'DPwi-'^)^'GDPw(^)) 

71 

{Ol_sp{x)Os-sFm 

^ [E(<' - ^n)' + {d( + dir] F^ix)] (A9) 



s-SF 

OC exp 



(O^^^^^(a:)O^i^y(0))cxexp 
(O^plW(^)O^P-5^^(0))«exp 



- 2 



F,{x) 



(AlO) 

'2[Yidif]F^{x)_, 

n 

(All) 



As a result, the correlation of different order will decay 
as a power law of x: 



(^G^''pwi^)^'GDPw(^)) 


oc 


x~ 


-2+aGDPW 


{oWwi^)o'c'Dwm 


OC 


x' 


-2+acDW 


{0'r'-SF{^)0','-SFm 


OC 


x~ 


-2+ap-sF 


{0's%Fi^)0's'^SFm 


OC 


x~ 





(A14) 



where 



CfGDPW 
CUCDW 
CUp^SF 

as-SF 



n 

2E(^^)': 



s 

^[J2{4 -4)' + id{ +dir]A15) 



Note that for CDW and p-SF phases, we just need to 
calculate the correlation between order parameters in 
the same tubes, because there is no correlation between 
tube j and tube j' in the definition of order parameters. 
Therefore, after we use the generalized Bogoliubov trans- 
formation to diagonalized the bosonic Hamiltonian of the 
Luttinger model, Eq. (|13p . the obtained transform ma- 
trix (Eqs. (jA4[) and (|A5p ) can be used to calculate the 
Luttinger exponents directly as shown above. 
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